Structure and Stability of Keplerian MHD Jets 



Thibaut Lery 1,2 , Adam Frank 3 
1 Department of Physics, Queen's University, Kingston, ON K7L 3N6, Canada 
2 Dublin Institute of Advanced Studies, 5 Merrion Square, Dublin 2, Ireland 
3 Dept. of Physics and Astronomy, Univ. of Rochester, Rochester, NY 14627-0171 

ABSTRACT 

MHD jet equilibria that depend on source properties are obtained using a 
simplified model for stationary, axisymmetric and rotating magnetized outflows 
(Lery et al. [Lery et al. 1998| , [Lery et al. 1999a] ) . The present rotation laws 



are more complex than previously considered and include a Keplerian disc. 
The ensuing jets have a dense, current- carrying central core surrounded by 
an outer collar with a return current. The intermediate part of the jet is 
almost current-free and is magnetically dominated. Most of the momentum is 
located around the axis in the dense core and this region is likely to dominate 
the dynamics of the jet. We address the linear stability and the non-linear 
development of instabilities for our models using both analytical and 2.5-D 
numerical simulation's. The instabilities seen in the simulations develop with 
a wavelength and growth time that are well matched by the stability analysis. 
The modes explored in this work may provide a natural explanation for knots 
observed in astrophysical jets. 



Subject headings: ISM:jets and outflows - MHD - Stars: formation - stability 



1. Introduction 

Observations of collimated outflows from young stellar objects (YSOs) and active 
galactic nuclei (AGN) suggest that magnetic fields may play a central role in the physics of 
these phenomena. Despite the large database of observations, theoretical MHD approaches 
have not yet converged on answers to several fundamental questions regarding the 
acceleration, collimation and propagation of these jets (See Lery et al. [Lery et al. 1999a| and 



reference therein). Difficulties investigating the nature of the jets arise mainly due to the 
high level of nonlinearity in the governing MHD equations particularly at critical points. 
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The work presented here is based on a simple model for jet launching and collimation 
presented in previous articles (Lery et al. [Lery et al. 1998|jLery et al. 1999a| , hereafter Paper 



I and II). It yields asymptotic MHD jet equilibria that account for the properties of the 
emitting source. The model is axisymmetric and stationary and in the work presented 
here includes more multi-component rotation laws for the source. The model assumes the 
magnetic surfaces possess a shape which is known a priori inside the fast critical surface, 
(we refer to the space within the fast surface as the inner region) As a first approximation 
magnetic surfaces were taken to be cones. 

The general problem of determining the stationary two-dimensional structure 
of magnetohydrodynamic outflows requires the solution of the equilibrium of forces 
perpendicular and parallel to the magnetic surfaces. One can describe the former by using 
the transfield or Grad-Shafranov equation, and the later using the Bernoulli equation 
for a polytropic equation of state. Asymptotically, the jet is assumed to be in pressure 
equilibrium with an external medium whose properties are independent of distance. The 
transfield and Bernoulli equations are solved both in the inner region and in the asymptotic 
region, which we also will refer to as the cylindrically collimated regime. 

The stability of the resulting cylindrical equilibria is of major importance, firstly 
because instabilities could explain observational features of jets such as wiggles, knots or 
helical filaments, and secondly, because they could be globally disruptive for the outflow. 
Thus, in this paper, we investigate both the linear stability and the non-linear evolution of 
our MHD equilibrium jets. This study will provide a diagnostic for features that should 
develop in fully time-dependent simulations of jet propagation. Such simulations will be 
the subject of another paper ( [Frank et al. 1999|) . 

Since jet behavior may be directly related to the properties of the emitting source, the 
models presented here can provide a basis for a better understanding of jet interactions with 
an ambient medium as well as propagation and stability issues in the context of the nature 
of the source. The main goals of the present paper are to produce, analyze and study more 
realistic MHD outflows than simple "top-hat" beams (constant density and velocity) used 
in previous studies ( [Frank et al. 199§ , kJcrqueira et al. 1999| . pardiner et al. 1999| , [Stone fc 



Hardee 1999|) . . These new models can then be used as input for numerical simulations. 
How outflows behave on critical surfaces and in the cylindrical regime constitutes the 
subject of the first part of the paper. In §2, we briefly recall the main features of the model. 
Jet equilibria and their properties are presented in §3. Stability analysis and non linear 
developments of the instabilities for the corresponding equilibria are investigated using ideal 
MHD computations and numerical simulations respectively in §4. Finally, we summarize 
the results in 55. 



- 3- 



2. The Model 

The Given Geometry (GG) model is based on the assumption that magnetic surfaces 
possess a shape which is known a priori inside the fast critical surface. For full details 
on the model and on justifications of the approximations, we suggest the reader see see 
Heyvaerts fc Norman 19891 , and Paper I and II. Inside the fast critical point the three 



critical points are assumed to be aligned on conical magnetic surfaces. Moreover the 
angular distribution of magnetic flux is assumed to be uniform. These approximations will 
be relaxed in forthcoming works. 

Unlike the Weber-Davis ( |1967| ) type models, however, the balance of forces 
perpendicular to magnetic surfaces is accounted for on the Alfven surface and at the base 
of the flow. On the Alfven surface the force balance equation becomes what is known as 
the Alfven regularity (non-singularity) condition. Owing to axial symmetry, stationarity, as 
well as the flux freezing condition, there exist in the general case five integrals of motion 
that are preserved on any axisymmetric magnetic surface (denoted a). Two of the integrals, 
the angular velocity Q(a) and a factor related to the entropy Q(a), are given as boundary 
conditions in the model. These two input fuctions constitute assumptions about the nature 
of the source rotator. 

MHD flows have two other critical points which are brought about by the Bernoulli 
equation. There are slow and fast magneto-sonic points both of which are located where 
the poloidal velocity equals one of the two magneto-sonic mode speeds. The fast and slow 
surfaces, unlike the Alfven surface, are saddle points, i.e., transonic solutions only exist 
for a certain relation among the integrals of motion. These relations are obtained from 
the criticality conditions that correspond to the vanishing of the differential form of the 
Bernoulli equation at constant a with respect to p and r. 

The Alfven regularity condition together with the criticality conditions determine 
the three other unknown integrals: namely the specific energy E(a); the specific angular 
momentum L(a); the mass to magnetic flux ratio a(a). We note that once these integrals are 
determined the model gives only an approximate solution because the transfield equation is 
not solved everywhere, but only at a few special places in the inner region. 

We further assume the density p to be related to the pressure p by a polytropic 
equation of state, p = Q{a)p 1 where 7 is the polytropic index. This assumption replaces 
consideration of energy balance and is meant to represent simply some more complex 
heating and cooling processes (See, for example, Vlahakis & Tsinganos 1998| for more 
general equations of state). 



Since the flow eventually becomes fully collimated far from the source, we use 
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cylindrical coordinates (r,(f>,z) and assume axisymmetry. R will be the spherical distance 
centered on the wind source. Each flux surface is labelled by the flux function a(r, z) 
proportional to the magnetic flux through a circle centered on the axis passing at point r, 
z. The physical flux is 2na. We normalize the magnetic flux to the total magnetic flux A 
enclosed in the jet, so that a is set to unity on the outer edge. It is convenient to split the 
magnetic field and the velocity into a poloidal part, which is in the meridional (r,z) plane, 
and a toroidal part. The former is denoted by a subscript P while the latter is just the 
azimuthal component. We note that the outflow does not necessarily fill all space from pole 
to equator. 

In the inner region, the system is governed by seven equations: the Alfven regularity 
condition; the four criticality conditions defined at the slow and fast magneto-sonic points; 
the Bernoulli equation written at the fast and slow critical points. We note that the 
number of equations can be reduced to six, due to the conservation of the energy on 
magnetic field lines. Using these equations the acceleration of the flow is studied up to the 
fast critical surface. When the flow exits the inner region, the flow begins to collimate 
as field lines bend towards the axis. We do not explicitly follow the collimation process 
between the fast critical surface {inner region) and the cylindrically collimated region. We 
note that between these regions a redistribution of the energy and angular momentum 
might occur due to shocks or other dissipative process. Nevertheless, the integrals of 
motion are assumed here to remain exactly the same all along the outflow for simplicity. 
Therefore, within the framework of the Given Geometry model, the asymptotic flow is 
uniquely determined from first-integrals obtained from the sub-fast surface regions. In the 
cylindrically collimated regime, the system consists of the asymptotic forms of the Bernoulli 
and transfield equations. In this region, z goes to oo, the flow density p is smaller than the 
Alfvenic density p\, r can be considered to be larger than r^, and the gravity becomes 
negligible. The full system of equations remains the same as in papers I and II and is given 
in Appendix A. 

We note again that in the inner region, the problem is entirely specified by two 
functions of a, namely Q(a) and Q(a), and also by one constant, the mass-to-magnetic flux 
ratio on the axis otQ. In the asymptotic regime, only one parameter is needed which can be 
either the external confining pressure p ex t, or the axial density po- 

Given the properties of the central emitting object, the model allows one to compute 
the dimensionless rotation parameter uo = Qr^/vpA, where r\ = L/Q is the Alfven 
radius and vpa the Alfvenic poloidal speed at the Alfven point vpa = (a|Va|/ (pr)) A . 
The density at the Alfven point is given by pa = Poet 2 and the Alfvenic Mach number 
by M\ = vp/vp A = pa/ P- More generally the subscript A will refer to values at the 
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Alfven point. Quantities referring to the two other critical points, namely slow and fast 
magneto-sonic critical points, will be indicated by subscripts s and f respectively. Rotators 
can be defined as slow, intermediate or fast according to whether the maximum of uo is 
much less than unity, close to unity, or near its maximum value for of (|) 3 ^ 2 (See Paper II 
for details). 



2.1. The Rotation Laws 



It is now well known that outflows can be accelerated from an accretion disk ("Disk 
wind" , Blandford & Payne gggg , Pelletier & Pudritz Fiege & Henriksen [Tggg , 

Contopoulos & Lovelace |1994j , Ferreira & Pelletier |1993| , Ferreira |1997| , Vlahakis & 
Tsinganos |1998| , Lery et al. |Lery et al. 1999b|) , or at the disk-magnetosphere boundary 
("X-winds", Shu et al. |1988| , |1994j ). Therefore, a Keplerian rotation profile is one of the 
most realistic description for rotation of sources producing jets. In paper I and II, we 
have assumed a rigid body rotation. Here, we have chosen more realistic rotation laws 
by including the Keplerian rotation in the outermost part of the outflow. In Fig. [I], the 
different profiles of angular velocity used in our study are presented. The pure Keplerian 
rotation law (dashed line) starts with a constant rotation close to the axis, as in the rigid 
body case (dot-dashed lines), but then follows a Keplerian profile. The multi-component 
(solid line), or multi-fast, case also starts with a rigid rotation corresponding, for example, 
to an axial ordinary wind. The angular velocity then doubles its value in order to model a 
source rotating more rapidly than the star in an intermediate region between the ordinary 
wind and the Keplerian disc wind that follows. Note that the angular velocity is always 
sub- Keplerian in the intermediate region. For all the rotation laws, the axial value of the 
angular velocity Qq is set to unity in the figure, and the radius is normalized to the size of 
the jet. In order to compare fast and slow rotators, we will respectively refer to multi-fast 
and multi-slow cases, the latter one rotating four times more slowly. 



3. The Jet Equilibrium 

3.1. The Numerical Procedure 

In the inner part of the flow, the variables calculated in the numerical procedure are 
the radii r s , r/, r^, and densities p s , pf, pa at the three critical surfaces along with the 
total energy E. In its cylindrical regime, the jet is entirely defined by r and p. All the other 
physical quantities can be derived from this set. For the numerical calculations, equations 
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have been reformulated as ordinary differential equations or converted from algebraic 
conditions into ODEs as functions of the flux surfaces a. 

The system consists of eight differential equations, namely the four equations of 
regularity on slow and fast surfaces (two on each surface), the Alfven regularity, the 
conservation of energy (a total of six for the inner region), and the Bernoulli and transfield 
equations in their asymptotic form (for the cylindrically collimated regime). 

The numerical solutions are obtained by initiating the integration of the system from 
the axis. Given the input parameters Qq, Oq, «o, 7, and po, all the critical positions 
and densities can be numerically obtained using analytical formulae (see Paper I). For 
numerical convenience, we prefer to provide the axial asymptotic density po rather 
than the external pressure. We further constrain the solution to be super-Alfvenic and 
super-fast-magnetosonic on the axis in the asymptotic region. This gives a limiting range 
of variations for input parameters and particularly to the axial density po. 



3.2. The Input Parameters 

The input parameters of the model can be selected so as to qualitatively reproduce 
observations. Given the properties of the jet-emitting object, i.e., its radius i?*, its 
temperature T*, the total mass loss rate M*, the base density n*, the magnetic field -B*, the 
factor and 7, it is possible to deduce dimensionless parameters f2, Q, olq. The parameter 
a>a can be a-posteriori related to the mass loss rate M#, i?*, and t he magnet ic field B*. Thus 
we define = ZkT^n*/ {mpn*)" 1 , a* = M^/AixR^B*, and Q* = ^GM^/Rf. All quantities 
are non-dimensionalized to reference values by setting Q = Q*/Q re f, «o = (x*/at re f and 
Q = Q*/Q re f. The entropy factor Q(a) is assumed to be constant across the jet. We have 
also studied models where the jet is hotter along the symmetry axis than on its outer edges 
and find that the results do not significantly change. 

In the present paper, we have chosen to model YSO jets with different rotation laws 
using typical values for T-Tauri stars as presented by Bertout et al. ( |1988| ) with M* = O.8M 
and = 3R Q . At the base of flow, we deduce the corresponding dimensionless input 
parameters: Q = 0.87, Q = 2, «o = 0.7, and po = 5.10 -7 . Major quantities of reference 
are then given (in CGS) by R re f = 10 15 cm, p re f = 250 ppc, v re f = 10 7 cm s -1 for Young 
Stellar Objects. The model could be applied to relativistic galactic sources or to quasars, 
by adapting the reference values and replacing the external pressure of the medium by the 
inertia of an electro-magnetic field outside the light cylinder in the relativistic regime. 
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3.3. The Critical Surfaces 

The system is solved on the three critical surfaces in the inner region. It has been 
shown in Paper I that the properties of these surfaces are related to the rotation parameter 
uj. As the parameter increases, the fast and Alfvenic surfaces move apart. In the lower 
right panel of Fig. |], the rotation parameter is plotted as a function of magnetic flux. 
The corresponding input parameters have been given in the previous section. This figure 
shows that the rotation parameter vanishes at the boundaries and reaches a maximum of 
1.6 around a = 0.45. Thus the jet generated by this flow will have the properties of a fast 
magnetic rotator at the location of the maximum in uj. Another relevant quantity aE/Q is 
plotted with the rotation parameter. Its minimum coincides with the maximum in uj. One 
can show that if this quantity possesses a minimum below its value on the polar axis, the 
flow will collimate cylindrically (at least for the region inside the minimum). This result 
is in good agreement with Bogovalov and Tsinganos (|1999[) who have shown that there 



always exist field lines that will collimate cylindrically for rotating MHD jets. In our case 
the cylindrical collimation is ensured due to the pressure of the external medium. 

The three critical surfaces and their corresponding densities are represented on Fig. 0, 
in the upper and lower left panels respectively. The surfaces are distorted with the largest 
deformation occurring where the rotation parameter reaches a maximum. A relation 
between the fast and Alfvenic radii has been found in paper I for the fast rotator case that 
states 
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This expression shows that when uj approaches its maximum value or when (3 vanishes, the 
fast point is pushed far from the Alfven point. It can even be rejected to infinity in the 
limiting cases. This is a well known result in the cold plasma limit (Kennel et al. ( |1989| )). 
In our case, f3 reaches 10~ 2 and u = 1.6 at the point where the critical surfaces are the 
most elongated. This behavior has also been obtained in Sakurai (|1987|) and Belcher & 
McGregor ( p76|) . 

When rotational effects are small (u ~ 0), i.e., close to the axis, the previous expression 
becomes r^/r/ oc 1 — 2uj 2 (slow rotator case in paper I). Thus the fast point gets closer to 
the Alfvenic point in these regions. 

The slow surface has the opposite behavior compared to the fast and Alfvenic 
surfaces. It deflates and gets closer to the source as uj increases, as one would expect from 
r s /rA oc (r^/77) 1 / 3 (Paper I). In fact in the fast rotator regime, the slow mode speed gets 
closer to the sound speed, and the slow mode acquires the character of a sound wave guided 
along the field line. Thus the fast magneto-sonic and Alfven surfaces strongly inflate when 
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rotation increases or the flow becoming cold. On the other hand, the slow surface gets 
smaller with increasing u>. 

The lower left part of the figure clearly shows the trend of decreasing densities from 
the axis to the maximum in u. At the fast point for fast rotators (see paper I), 

EL ku -^(1A)\ (2) 

PA Vf) 

The square of the last term explains the decrease of four decades in density. As stated 
earlier, all the other variables of the problem, such as the first integrals of the motion E, L 
and a can be deduced from the positions and densities of the critical points. For example, 
Fig. |2| shows that the total angular momentum reaches a maximum inside the jet and not on 
the outer edge. Thus most of the dynamics will be internal since this quantity is conserved 
along the flow. 



3.4. The Asymptotic Equilibrium 

3.4-1- The Numerical Solutions 

Here we compare the influence of the rotation laws on the cross-sectional distributions 
in MHD jets. The quantities that define the jet in the cylindrically collimated regime 
are plotted in Fig. |3| for pure Keplerian, multi-fast and constant rotations. The z and (f) 
components of velocity and magnetic field are represented together with the density p and 
the net electric current la, as functions of the relative radius (normalized to the jet radius). 
The length scale is the jet radius. The density is normalized to its value on the jet axis Pq. 
The non-dimensional velocities refer to the fast magnetosonic velocity = <? s + v\ on the 
axis, c s being the sound speed. The magnetic field is normalized to v /po Vf. 

The poloidal velocity v z increases from the axis to the outer edge for the rigid rotation 
case, and slightly decreases for the pure Keplerian case. For the multi-fast case, however, v z 
peaks where the rotation parameter is maximum, approximately in the middle of the jet. 
Therefore in the latter case, the fastest part of the jet is neither on the axis or on the outer 
edge, but inside the jet itself. The azimuthal velocity follows the same trend but with 
a magnitude several orders smaller than the poloidal component. If v z is of the order of 
several hundreds of km s -1 , is only about a few km s _1 . 

The azimuthal component of the magnetic field always dominates the poloidal part 
except at the axis. Thus for these parameters the jet's field is highly twisted. The azimuthal 
component follows approximatively a 1/r law in the outer part. 
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As one moves towards the axis from the edge of the jet magnetic pinch forces or hoop 
stress become increasingly important. In order to maintain an equilibrium the gas pressure 
must balance the hoop stresses. Thus we see large pressure and density gradients in near 
the axis. We denote the high density region centered on the axis as the core, and the lower 
density outer regions as the collar. Note that the bulk of the jet's momentum resides in the 
core. Hence this portion of the beam will penetrate more easily into the ambient medium 
during the jet's propagation while the collar will be more strongly decelerated. Thus we 
expect that even if the relative velocity is smaller close to the axis, the central part of the 
flow will propagate faster. 

The last quantity represented in Fig. ^| is the "poloidal electric current" through a 
circle centered on the axis and extending out to a magnetic surface a. This quantity is given 
by Ic(a) = —rB^j \iq (The physical poloidal current is I p h ys = — 2irlc). While the current is 
always increasing for a constant rotation jet, in the Keplerian and multi-component models 
it reaches a maximum in the middle of the beam and vanishes at the outer boundary. 
Therefore a return current flows back inside the jet in these cases. This will be studied in 
more detail for the multi-fast case in the next section. 

The kinetic to Poynting flux ratio (not shown here) can be calculated as a function 
of the relative radius. We find that only the central part of the asymptotic outflow has 
a kinetic energy flux strongly dominating the Poynting flux. Thus away from the axis a 
non-negligible part of the magnetic energy is not transformed into kinetic energy. 

Thus rigid body rotation jets are characterized by a dense, current- carrying core 
having most of the momentum, surrounded by a tenuous current-free envelope, dominated 
by the azimuthal magnetic field. Keplerian and multi-fast cases also exhibit a central 
current- carrying core but are surrounded by a denser collar and carry an internal return 
current. 



3.4-2. An Approximate Analytical Solution 

By simplifying the asymptotic equations, it is possible to obtain approximate analytical 
solutions in the region where the jet is cylindrical. Away from the axis, where the gas 
pressure is negligible with respect to magnetic pressure, the transfield equation can be 
integrated to yield f2 2 r 4 p 2 /iQ or 2 = constant. Then the density can be expressed as a 
function of r and of the first integrals 



p(r) « Ca/nr 2 , 



(3) 
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where C is a constant. In the present work, the angular velocity Q(r) is given as an initial 
condition, while the mass to magnetic flux ratio a(r) is part of the solution. If a(r) and 
Q(r) are constant, the density drops as r~ 2 . This is the behavior that we get for the rigid 
body rotator away from the axis. On the other hand, for increasing a, or decreasing values 
of Q, the density can increase as obtained for the multi-fast rotator in the outer regions 
of the jet. Close to the axis, i.e., r <C 1, the gas pressure dominates and the density 
approaches a constant, the transfield equation reducing approximately to Qp 1 = constant. 

The above result allows one to derive similar approximate formulae for the toroidal 
components of the magnetic field and velocity which respectively reduce to \B^,\ ~ Qpr/a, 
and \v,p\ ~ pflr / \IqO? . Combined with Eq. ^, the equations become 

|2J,(r)|«Cyr , Mr)\*C/a. (4) 

It follows that the net asymptotic electric current is approximately given by 
/ = 2-nr\B^\j 'p « constant. Therefore, in the region where magnetic pressure 
dominates over gas pressure, the asymptotic current will be constant. Finally, one can 
deduce the asymptotic poloidal velocity of the flow v z (r) = (a/ pr)da/dr, that becomes 

v z (r) « fir/C. (5) 

Fig. |] shows the approximate analytical solutions together with the solution for 
the multi-fast case. The density and the poloidal velocity are well reproduced in the 
intermediate region where the magnetic pressure dominates. The toroidal component of the 
magnetic field is well described by the formula only far from the axis. Finally, the toroidal 
component of the velocity is similar in the analytical and the numerical solutions. 

In theory, by using this approximate solution, it may be possible to deduce the 
properties of an outflow and, therefore, of the emitting object directly from quantities 
observed in the jets. A rough estimate of a could be deduced from the toroidal component 
of the velocity. Combined with an accurate measure of the density and of the asymptotic 
axial velocity, it would give the angular velocity. Using the model, the five first integrals 
could be estimated and would allow a characterization of the source itself. Moreover, just 
from the profile of the density, it would be possible to deduce the mass to magnetic flux 
ratio, a(r), for a given rotation law. Thus a precise measure of the mass density in jets may 
be a key point in order to deduce the properties of the source at its base. 



- 11 - 



3.5. The Current along the Jet 

In Fig. [|, we present the net electric current on the Alfvenic and fast surfaces and in the 
cylindrically collimated region, in order to have a global picture of the current circulation in 
the jet. Only the multi-fast case is presented here since the Keplerian case presents similar 
trends. The current first increases outwards from the axis, then reaches a plateau where 
the magnetic pressure dominates, as shown in the previous section. Finally, it decreases in 
the outer part of the jet, and almost vanishes on the outer boundary. The current density 
is negative in the outer part of the jet and therefore the direction of the current is opposite 
to the one in the axial region. Thus there exists a strong current in the core and a return 
current in the collar, the intermediate part of the jet being almost current-free. 

We note that the current is never strictly zero on the outer boundary, and therefore 
a small part of the current may flow back either on the external surface of the jet, or 
in the external surrounding medium. This result is similar to the force-free field model 



of Lynden-Bell (|1996|) , where the current circulates inside the jet itself, and where the 



magnetic field lines are anchored in a differentially rotating accretion disc. 



3.6. Comparisons with Previous Works 

The characteristics of magnetized outflows at large distances from the central object 
have been addressed by a variety of studies. Many results of the present study are in good 
agreement with these works. For example, as found by Beskin et al. ( |1998|) , a jet with a 



zero total electric current has its angular velocity vanishing at the jet boundary, and its 
axial regions dominated by kinetic energy. Concerning the current, Bogovalov ( |1995| ) has 



shown analytically that there always exists a field line in the outflow which encloses a finite 
total current, and consequently, the asymptotics always contain a cylindrically collimated 
core. This is precisely what we find in the axial region of our jets, regardless of the presence 
of any external medium that would ensure cylindrical collimation anyway. As with Ferreira 
( |1997|) , it is possible to show, as it is done in paper I, that the minimum mass loss rate has 



a lower limit and can not be arbitrarily small. We also agree with Ostriker ( Pstriker 1997 ) 



and Lery et al. ( |Lery et al. 1999fj| ) who conclude that the optical jet may represent only the 
densest part of the total outflow. 

We obtain fast magnetosonic Mach numbers, (which also corresponds to the Alfvenic 
Mach number on the axis), between 2 and 4. This range corresponds to what Camenzind 
( |1997|) has found for his model for low-mass protostellar objects. The corresponding 



jets have low fast magnetosonic Mach-numbers Ma — 2. By taking into account an 
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accretion disc around the stellar magnetosphere, Fendt & Camenzind ( |1996| ) also find a fast 
magnetosonic Mach-number m 2.5. However, this does not seem to be a general statement 
about MHD jets since there exist models with larger values (Sauty et al. ( |1994| ), Trussoni 
et al. (TO)). 

Shu et al. ( |1995| ) have studied magneto-centrifugally driven flows from young stars 
including their structure at large distances from the source when they collimate into jets. 
When R ^> 1 and r is at least moderately large also, they have found that the density 
distribution can be given by p — > K(r)/a(r)r 2 (their Eq.4a) where K(r) is an arbitrary (but 
slowly varying) positive function of r. This result, and the variations of their asymptotic 
velocity, are in good agreement with Eq. ^, with K(r) being C/Q(r) in the present model. 
Despite the similarity of the analytical results, the functions in Shu et al. do not correspond 
to the multi-fast (nor the Keplerian) case. Note that our model does not describe the 
physical processes occurring at the source itself, i.e., at the surface of the disk or the 
disk-star boundary. Thus one may be able to link our results to models such as Shu's for 
the generation of the outflow at its source. 



3.7. The Astrophysical Consequences 

The model makes it possible to obtain MHD jet equilibria with multiple components, 
i.e., a central dense core surrounded by a collar. Characteristic physical quantities are close 
to those obtained observationally, with, for example, velocities of several 10 2 km s^ 1 for 
YSOs. We note that there exists a broad range of solutions to the equations that possess 
favorable characteristics for compararisons with observations. Moreover, the properties 
of the different parts of the jet can be directly related to the properties of the source. 
Therefore the ensuing equilibria can be used to simulate numerically the propagation of 
the jet into an external medium, more realistically than with a "top-hat" distributions for 
velocity and density. By comparing astrophysical observations with the results of different 
simulations of jet propagation could ultimately provide a tool for deriving source properties. 

Another important issue that must be addressed is the effect of ambipolar diffusion in 
our model, and the estimation of the distance where it should become effective. In a recent 
work by Frank et al. ( |1999| ), it has been shown that the ambipolar diffusion time-scale was 



related to the plasma ft parameter, ft = P 9 /Pb, by 



T0 3 cm- 3 / V10 15 cm/ \ T jet J \ft + l 
In this equation, the ambipolar diffusion time-scale is given in years, n n is the number 
density of neutral particles, and Tj et is the temperature. When ft — > oo, the term ft /(ft + 1 
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tends to unity, while when < 1, the ambipolar diffusion time-scale can be effectively 
reduced. In the present model, the plasma (3 parameter is given by 

0(a) = 27 9Ea1 ( 7 ) 
f [ ) (7 - 1) W ■ [7) 

As shown in Fig. |6], (3 can become quite low in some part of the jet. The most dramatic 
case is the multi-fast case, where value around 10 -2 can be reached in the intermediate part 
of the jet, the core and the collar having a plasma parameter larger than unity. Therefore, 
the jet is magnetically dominated only in the intermediate zone, precisely where ambipolar 
diffusion can be effective. For all the other cases, the average value of the parameter is 
about unity. From observations, typical values for YSO jets are: 10 3 < n n / cm~ 3 < 10 4 ; 
1 x 10 15 < rjetl cm < 5 x 10 15 ; 5 x 10 3 K < T < 3 x 10 4 K (Baccioti & Eisloeffel p99l) . 
For jet parameters in the middle of the expected range of variation, we find t a d of order 
10 5 to 10 4 y. Thus the dynamical time-scale for YSO jets is t dyn = 10 4 — 10 5 y (Reipurth 
et al. |1997| , Eisloeffel & Mundt |1997| ). This is of order of, or greater than, the ambipolar 
diffusion time: td yn > t a d- If we consider length scales, the distance D a d, where ambipolar 
diffusion becomes effective is approximatively between 0.1 pc and 3 pc, for jet velocities 
around 10 2 km s . This is the range of distances where the magnetic field might be altered 
or at least where internal configurations should evolve. Such an effect should reduce the 
internal hoop stress that would normally help the jet to remain collimated. It could also 
reduce instabilities as it will be shown in the next section. Hence, parsec scales YSO jets 
should naturally be less collimated when td yn ~ 10 4 — 10 5 y. 

Finally, the present jet equilibria show large gradients in velocity, magnetic field and 
density. Such configurations should give rise to instabilities that could change internal 
structures of jets or even disrupt them. This is the subject of the next section. 



4. The Stability Analysis 



Current-carrying jets with a helical fields are susceptible to pressure-driven (PD), 
Kelvin-Helmholtz (KH) and magnetic instabilities driven by the electrical current (Current 
Driven, CD). The latter modes are due to field aligned electrical currents. Generally 
instabilities will be a mixture of CD, PD and KH modes with the contributions from each 
type of mode being difficult to estimate. Stability properties of such MHD jets have been 
investigated only recently (Appl & Camenzind |1993| , Appl |1996| , Lery [Lery 1996| , Appl et 



al. |1999| ). These instabilities could play an important role in various observed morphological 
structures such as wiggles, (e.g. for quasars: Krichbaum et al. |1990| , Feretti et al. |1999| , 
for YSO jets: Schwartz & Greene |1999] ), knots (e.g. for YSO jets: Ray et al. |1996| , Raga 
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& Noriega- Crespo |1998| , Rosado et al. |1999|) , and filaments {e.g. for AGNs: Biretta |199€ 



Bahcall et al. [I995| , for YSO jets: Elmegreen |1989| , Dutrey et al. |T99T 



4.1. The Linear Stability Analysis 

In order to study linear development of instabilities, we adopt the standard (temporal) 
approach where the axial wavenumber k is real and the imaginary part of the complex 
frequency corresponds to growth rate T, growth time being r = Wavenumbers are 
given in units of inverse jet radius and growth rate is normalized to the inverse Alfven time. 
A synopsis of the method is presented in Appendix B. 



4-1.1. The Pinch Mode m = 

Stability analysis of magnetic configurations derived from the Given Geometry model 
have been studied in Lery ([Lery 19^5]) for rigid rotators. We report here the main results. 



It has been shown, in such a case, that the pinch, or sausage, mode (m = 0) dominates over 
other m-th order modes. This is mainly due to the large gradients of the density and of the 
magnetic field that create strong PD and CD instabilities respectively. 



The results of (|Lery 1996| ) show that pure CD modes can develop on rapid, i.e., Alfven, 



time scales. Moreover fast rigid rotators are more unstable than slow ones. In Fig. [7], we 
have plotted the dispersion relation for these two cases for the pinch mode m = 0. By 
studying a large range of angular velocities, Lery (|Lery 1996|) demonstrated that faster 



rotation rates produce larger growth rates and smaller wavenumber cut-offs. This is 
clearly seen in the figure, where the fast rotator is the most unstable for large wavelengths 
(i.e., small wavenumbers), but also presents a smaller cut-off in wavenumbers. The 
maximum growth rates and their corresponding wavenumbers are reported in table [l|. The 
unstable modes that should grow the most rapidly have a wavelength about 3 jet radii for 
the fast rotators and half the jet radius for the slow rotators. These results can be applied 
to observations of real YSO jets. For HH34 the largest wavelength above corresponds 
approximately to the minimum knot separation of 3.4 rj et as given by Burke et al. ( 198£| ) . 



In the case of HH111 however, the separation is larger, about 11 Tj et (Reipurth |1989| , Morse 



et al. |1998|) . Thus model seems not to fit to the HH111 case. If the jet diameter is not 
defined as the full width of the jet, 2 rj et as we have done, but as the observationally 
deduced knot width then the application appears better. This would correspond, in our 
model, to the denser, axial core. In such a case, the jet diameter would be smaller by a 
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factor of 3 or 4 and knot separation would be about 9 to 12 r jef , thus accounting for the 
larger separation of knots of HH111. 

We note also that the number of knots for HH111 is about 13 (Reipurth 
1989|) corresponding to an approximate length of the knotted section of the jet of 



4 x 10 17 cm = 0.13 pc, (given a knot separation of 3 rj et = 3 x 10 16 cm. This is the distance 
where ambipolar diffusion could become effective in the jet, as pointed out earlier in the 
present paper. Therefore, at 0.1 pc from the source, the internal configurations should 
begin to evolve leading to decreases in pressure and magnetic field gradients. Consequently 
pressure driven and current driven instabilities should be less important. This might be the 
reason for the disappearance of the knots at this distance from the source. 

Thus, regardless the type of rotators, the domination of axisymmetric modes and their 
corresponding wavelengths suggest that these instabilities could be at the origin of the knotted 
structures seen in a large number of jets as seen, for example, in HL Tau, HH1, HH30 and 
HH34 (Ray et al. \IUUflj) , or in HH11 and HH311 (Rosado et al. \199$ ). 



This potential origin for knots has already been proposed by several authors (Burke et 
al. |1988| , Bodo et al. |1995| , Micono et al. |1998|) . These studies focused on Kelvin- Helmholtz 



instabilities. Todo et al. (|1992|) , however, have argued that the pure KH instabilities could 
not alone produce the knots. We propose that the knots may be due to the combined effects 
of the KH, CD, and PD instabilities, with a domination of the PD and MHD instabilities. 



4.1.2. The Helical Mode \m\ = 1 

As stated above an MHD outflow can exhibit Kelvin-Helmholtz, current-driven, and 
pressure- driven instabilities. In general, it is not possible to entirely separate these modess. 
Only with very simple and idealized configurations one can hope to understand the physics 
driving individual modes. In this section we restrict ourselves to current- driven modes, in 
order to study their effects on our models. We use the results obtained by Appl et al. ( |1999| ) . 
They have approximated the MHD jet by an infinitely long cylindrical outflow of a perfectly 
conducting plasma. The jet is supposed to have constant density and velocity, as well as 
negligible thermal pressure and rotation. Consequently, Kelvin-Helmholtz instabilities arise 
only due to the vortex sheet at the jet boundary. Moreover, pressure-driven instabilities 
are excluded by considering cold jets. The simplifications differ from our present model, 
i.e., we have velocity and pressure gradients. Nevertheless it allows to get a rough estimate 
of the maximum growth rate and wavenumber for each value of the axial pitch. Eventually 
the ensuing characteristic wavelengths can be calculated. 
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Appl et al. ( |1999| ) have shown that CD helical modes generally dominate in thier jet 
configurations. They found that the magnetic pitch, P = rB z j{r^ et B^)^ and in particular 
its value on the axis, Pq, essentially determines the growth rate of the helical instabilities in 
the case of small values of the pitch. We have plotted the pitch in Fig. [8] (left panel) as a 
function of the relative radius, for Keplerian, multi-slow and multi-fast rotators. The axial 
values of the pitch functions were deduced from the models presented above. Following 
the same procedure as Appl et al. (|1999|) , it is possible to derive the maximum growth 



rates and wavenumbers as functions of the axial pitch for our jets (right panel in Fig. [|). 
Results are reported in table |2|. The most unstable magnetic configurations correspond to 
the multi-slow case which has a very small wavelength. We note that the Keplerian and 
multi-fast cases have maximum growth rates that are still quite large. 

One should note that the maximum growth rates given here have been obtained for a 
particular value of medium external pressure. An increase of the external pressure increases 
the central pitch Pq which reduces T max by a factor of 10. For large external pressure, the 
pinch mode can be of the same order or larger than the helical mode. A 3D simulation 
would be needed in order to study the development and the evolution of these instabilities 
in detail. 



4.2. The Non-linear Evolution 

4.2.1. The Kink Modes 

Non-linear computations of the development of the CD modes have been carried out by 
Lery, Baty & Appl ( |Lery et al. 1999c| ) using a cylindrical evolution code. These calculations 



show that like the linear evolution, non-linear instabilities also develop on rapid time scales 
and therefore should affect internal jet structures. Moreover the non-linear behavior of the 
jet equilibrium is highly sensitive to the structure of the initial magnetic configuration. 
In the present case where the pitch function is increasing, the dominant mode (m = 1) is 
most probably resonant at half the relative radius of the jet. A current sheet should form 
at this locus producing a kinked inner cylinder due to helical distortion. Reconnection and 
turbulence in such a configuration should occur and accelerate particles on this surface. 
The ensuing saturated configurations may look like a hollow cylinder. For AGNs, such 
a structure may be related to the formation of the VLBI knots, which apparently are 
non-axisymmetric features within the flow channel (Krichbaum et al. |199U| ). Finally, a 
high current density forms along the jet axis in these calculations. This should give rise to 
radially localized dissipation which could be a potential heating mechanism within the jet 
core. 
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4-2.2. The Simulations 

The Numerical Method A multi-dimensional (2.5-D) simulation of a multi-fast MHD 
jet have been performed using a MHD TVD code in cylindrical symmetry (see Ryu, Jones 
& Frank 1995a and Ryu et al. 1995b|) . The simulations are initiated with the cylindrically 



collimated jet equilibrium traversing the length of the computational domain (256 x 1024 
zones), the jet radius being 64 zones. The equilibrium is then continuously injected at 
z = boundary of the grid. The radius of the jet was chosen to correspond to a YSO jet of 
10 16 cm. The maximum density of the initial jet was 250 ppc, and the central velocity was 
10 7 cm The pressure outside the jet was imposed by the model through the pressure 
balance at the outer boundary. The jet is surrounded, in the present simulation, by a 
magnetized medium that has a small poloidal field and no toroidal component. 

As seen previously, the present magnetic configuration is always unstable to both the 
pinch and the helical (or kink) modes. Since our simulation is axisymmetric, it is only 
possible to track the pinch modes. In order to study the kink modes 3D simulations are 
required. This will be the subject of a forthcoming paper. 



The Results We note first that the equilibrium remains stable over a long period. No 
preferential perturbation has been used in order to destabilize the simulation. As the 
jet propagates numerical noise generates the instabilities. Two types of regularly spaced 
features arise. One has a small wavelength and is located close to the axis. The second 
presents a larger wavelength, the second mode is wider and situated at the interface between 



the jet and the ambient medium as illustrated in fig. and fig. ITUI where the density, and 



the poloidal and toroidal components of the magnetic field are shown respectively. 

In a companion paper (Frank et al 1999) we have performed simulations which follow 
the propagation of the jet models presented here. These simulations also show instabilities 
developing in the beam which have identical wavelengths. The instabilities are initiated 
where the beam interacts strongly with ambient medium via the propagation. Here the 
instabilities take longer to develop as they are generated via low level noise. Note the form 
of the instability at the axis. Detailed inspection of the simulations show a strong axial 
flow (a pinch). When material reflects off the axis it expands in a tight balloon or loop-like 
structure. As the flow evolves each loop interacts with its neighbor leading to a saturation 
of the growth. Consideration of the poloidal and toroidal fields shows that initially the 
loops are dominated by the toroidal components. 

In Fig. |9], the gray-scale map of the density, we have drawn schematically the envelope 
of the largest instabilities (dashed line) and marked by large arrows the corresponding 
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location of the knots. The smallest internal instabilities are also shown by smaller grey 
arrows. We note that there exist globally five small instabilities between two consecutive 
knots of the envelope. The corresponding wavelengths are reported in table [| The 
corresponding wavenumbers are also provided. 

We note that the wavelength of the envelope is around 3r je4 and that the smallest 



instabilities correspond to half the jet radius. As shown by Appl et al. ( |1999| ), current 
driven instabilities are absolute instabilities, i.e., they grow but do not propagate, in the 
rest frame of the jet, and the boundary conditions at the outer edge of the jet do not affect 
the internal modes. Hence internal instabilities will develop first as seen in Fig. |9|. On 
the other hand, the instabilities with large wavelengths that develop on the envelope will 
directly depend on the the boundary of the jet that is only perturbed when the internal 
instabilities become important. Thus in real astrophysical jets, small internal instabilities 
should develop before larger ones which give rise to the larger pinches of the jet. 



4.3. Comparison of the Results 

Since our simulation is axisymmetric and therefore only track the pinch modes, the 
results of the simulation can only be compared with those presented in section |4.1.1| . 

In the axial part of the jet, the rotation parameter uj is small with respect to unity. 
Following the classification of paper I, this corresponds to the characteristic of a slow 
rotator. Hence, for our simulation, the instabilities in the core should be related to the 
results presented in table [I] for the slow rotator. On the other hand, in the intermediate and 
the outermost regions, the jet properties correspond to those of a fast rotator. The largest 
growth rate for the fast rotator corresponds to that of the envelope. These wavelengths 
have been reported in table f|, together with their ratio, and show good agreement. 

The first internal instabilities arise at around 5 rj et in the simulation. They grow and 
become important at about 10 r jet . The mean time for the instabilities to grow is then 
approximatively given by 

Tj « 5 r jet /v jet . (8) 

Let us compare this with the growth time of the pinch mode instabilities t s i ow for the slow 
rigid rotator. The inverse of the maximum growth rate for slow rotator (table [I]) gives it 
as Tsiow = Tj low = 1.5 tai where ta = rj et /vA is the Alfven time. Using the value of the 
relative velocity in Fig. [3] (upper left panel), we can derive the velocity of the jet w.r.t. the 
Alfven velocity, i.e., v = Vj et /vA ~ 3 ± 0.5. Then the Alfven time is approximately given by 
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ta ~ 3(±0. 5)rj et /vj et . Hence we finally have 

t s1ow = 1.5 t a « 4.5 (±0.75) r jet /v jet , (9) 

that can be compared to rj. Hence the instabilities develop in the simulation as predicted 
by the stability analysis. Moreover these growth times are rather short and consequently 
the internal instabilities in astrophysical jets should develop relatively close to the source 
(at about 5 rj et ). [] 

Thus the theoretical results are in good agreement with the simulations. The present 
stability analysis seems to be a good diagnostic tool to understand the development of 
instabilities in MHD jets, and to predict their characteristic wavelengths. 



5. Conclusions 

In this work we have investigated Keplerian MHD jet equilibria by means of a simplified 
axisymmetric, polytropic and stationary model. It assumes that the magnetic surfaces 
possess a shape, conical in the present case, which is known a priori inside the fast critical 
surface. In this inner region, the problem is entirely specified by Q(a), Q(a), and atQ. The 
Alfven regularity condition together with the criticality conditions determine the three 
other unknown first integrals: E(a), L(a) and a(a). In the asymptotic region, the flow 
is uniquely determined by the first-integrals obtained in the inner region together with 
the axial density po which is given as boundary condition. This region is governed by the 
asymptotic forms of the Bernoulli and transfield equations. Our principal conclusions can 
be summarized as follows: 

1. In the inner region, we find that the fast magneto-sonic and Alfven surfaces strongly 
inflate when rotation increases, while the slow surface gets smaller. The fast point 
remains at a finite distance for finite entropy flows, in contrast to cold flows. 

2. In the asymptotic region, jets have a dense and current-carrying central core where 
most of the momentum is located. Rigid body rotation jets are surrounded by a 
tenuous current-free envelope, dominated by the azimuthal magnetic field. Keplerian 
and multi-fast jets are surrounded by a denser collar and carry an internal return 
current. The intermediate part of the jet is almost current-free and is magnetically 
dominated. 



1 As shown by Begelman (1998), the distance R m in from the source where instabilities start growing in 
a non-relativistic, supersonic jet must be larger than Ma rj e t (in our case R m in > 2.5 rjet)- 
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3. An approximate analytical solution was derived which may make it possible to 
estimate the properties of the outflow and of the emitting object directly from the 
quantities observed in jets, e.g. the density, the components of the jet velocity. 

4. We have addressed the linear stability and the non-linear development of instabilities 
for the corresponding equilibria using both analytical and 2.5-D numerical simulations. 
Instabilities in the simulations develop with a wavelength and growth time that 
are well matched by the stability analysis. Regardless the type of rotators, the 
wavelengths of the axisymmetric modes suggest that these instabilities could be at 
the origin of the knotted morphology of a large number of astrophysical jets. 

5. The rather short growth times of the internal instabilities suggest that they should 
develop relatively close to the source, at about 5 r jet for YSO jets. 

6. We have shown that ambipolar diffusion should be effective on times scales of order 
10 4 to 10 5 y, or a length scale equal or larger than 0.1 pc, at least in some part of 
typical YSO jets. The consequent changes of the internal magnetic configurations 
might be the reason for the disappearance of knots at such a distance from the source, 
as observed. 

7. The relatively good agreement between the theoretical results and the simulation 
shows that the present stability analysis is a good diagnostic tool to understand the 
development of instabilities in MHD jets. 

We wish to thank Hubert Baty, Stefan Appl, Jean Heyvaerts for their help with this 
project. We thank Tom Gardiner and Guy Delamarter for useful discussions and Tom Jones 
and Dongsu Ryu for their help with the numerical code. This work was supported by NSF 
Grant AST-0978765 and by the University of Rochester's Laboratory for Laser Energetics. 
The project was also supported by an operating grant from NSERC of Canada. 



1. The so-called Alfven regularity condition is the particular form assumed by the force 
balance perpendicular to the magnetic field, or transfield equation, at the Alfven point. It 
can be written as 



A. The Differential Equations of the Model 
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where p is the slope of the solution of the Bernoulli equation at the Alfven point. The 
reader should see [Heyvaerts fc Norman 1989| , [Lery et al. 1998| and [Lery et al. 1999a] for a 
more detailed explanation of the the equations and methods 

2. The projection of the equation of motion on Bp yields by integration the Bernoulli 
equation, E(a) = v 2 /2 + (7/(7 — 1)) Qp 1 ^ 1 + G(r, z) — rQB e / p a. In the absence of 
MHD forces, this first integral would express the well known Bernoulli's theorem, i.e. the 
constancy of the sum of the kinetic, enthalpy and gravitational energy fluxes. The presence 
of the magnetic field introduces another energy flux, the Poynting flux, the fourth term in 
the equation. The total specific energy has a given value on a magnetic surface. This can 
be expressed in differential form as 



d_ 
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where the s and / subscripts refer to the slow and fast points respectively 



3. Every solution has to fulfill four criticality conditions defined at the slow and fast 
magneto-sonic points. These points are located where the differential with respect to p 
and r of the Bernoulli function B(r,p) = E(a) vanishes, i.e., rdB/dr = 0, pdB/dp = 0. It 
is convenient to convert the algebraic magnetosonic criticality conditions into differential 
equations. The critical points can be followed from one magnetic surface to the next by 
differentiating these criticality equations with respect to a: 
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4. The asymptotic form of the Bernoulli equation can be written as: 
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5. The asymptotic form of the transfield equation is given by: 
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B. The Stability Analysis Method 

A global normal mode stability analysis is used. The radial displacements, £ r , of the 
fluid elements from their equilibrium positions are of the form, 

£ r (x) = £r( r ) expz(m0 + kz — out) (Bl) 

and similarly for the remaining perturbed quantities. The linearized equations can be 
cast into a system of two first order ODE for Y = r£ r and the total perturbed pressure 
Z = Sptot = S(p + B 2 /8n). It is given by 

dY 

AS— = dY - rC 2 Z, (B2) 
dr 

and 

AS— = C 3 Y/r - C X Z, (B3) 

where A, S and Cj functions of the equilibrium quantities as well as the Fourier parameters 
u,k,m (Appl & Camenzind [L993|) . 

Jets propagate through a compressible medium of high conductivity. Since one is 
interested only in instabilities due to the jet itself, radially outgoing decaying waves for 
large radii define the boundary condition, 

P oc (Ar) (B4) 

for r — > oo. The are the Hankel functions of the first kind and A 2 = cu 2 /c 2 — k 2 , with c s 
the sound speed, if the ambient medium is assumed to be unmagnetized. The perturbed jet 
is assumed to remain in equilibrium with its surroundings. Regularity on the axis provides 
the other boundary condition, i.e., Y = for r = 0. 
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Relative radius 



Fig. 1. — Rotation laws for constant (dot-dashed), pure Keplerian (dashed), and multi- 
components, or multi-fast, (solid) models. Axial angular velocity is set to unity. 
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Fig. 2. — Multi-fast case : We show the critical surfaces in dimension-less quantities (upper 
left), and the densities on these surfaces (lower left). Three first integrals, namely the total 
angular momentum L (upper right), the rotation parameter u and the ratio aE/fl (lower 
right) are calculated on the Alfvenic surface as functions of the magnetic flux a. Input 
parameters are Q = 0.87, SI = 2 and = 0.7 (7 1). 
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Fig. 3. — Variations with the relative radius of the velocity V and magnetic field B 
components, of the density p and of the net electric current Ic in the cylindrically collimated 
regime. Constant (dotted lines), pure Keplerian (dashed), and multi-fast (solid) rotation 
laws are considered. 




Fig. 4. — Jet variable distributions for multi-fast (solid line) and approximate analytical 
(dashed line) solutions in the asymptotic region as functions of the relative radius. 
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Fig. 5. — Schematic representation of the variations of the net electric current along the jet 
(heavy solid lines), and its direction (arrows), at the Alfvenic (A), and the fast (F) surfaces, 
and in the cylindrically collimated regime (C). The light solid (and dashed) lines represent 
the poloidal projection of magnetic field lines. 




Fig. 6. — Plasma (5 parameter as a function of the relative radius for the multi-fast (solid), 
multi-slow (dashed), rigid body (long dashed) and pure Keplerian (dotted-dashed) rotation 
cases. 
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Fig. 7. — Dispersion relation of the pinch mode (m = 0) for fast and slow rotators (solid 
and dotted lines respectively). Dashed lines correspond to the maximum wavenumbers 
{kmax = 2.02, k max = 11.82). 
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Fig. 8. — The Pitch profile (left panel) is shown for the Keplerian (dotted lines), the multi- 
fast (solid), and the multi-slow (dashed) cases. On the right panel, the maximum growth 
rate (T max ) and the maximum wavenumber (k max ) are plotted as functions of the central 
pitch value. The arrows mark their estimate for the different cases. 



Fig. 9. — Structures of instabilities. Thick white arrows correspond to knots of the envelope 
(delimited by dashed line), and thin black arrows to maxima of density for the internal 
instabilities. 




Fig. 10. — Components of the magnetic field of the destabilized equilibrium 



-35 - 



Table 1. The Pinch Mode: growth rate T max , wave number k max , and wavelength \ max of 
the instability with the largest growth rate for the slow and the fast rigid body rotators. 



type 


r 

L max 


kmax 


^max 


Fast rotator 


3.50 


2.02 


3.11 


Slow rotator 


0.67 


11.82 


0.53 



Table 2. The Helical Mode: central pitch P , growth rate T max , wavenumber k max , and 
wavelength \ max of the instability with the largest growth rate for the different rotators. 



type 


Po 


r 

i max 


h 


^max 


Keplerian 


0.031 


4.2 


25 


0.25 


multi-fast 


0.023 


5.5 


30 


0.20 


multi-slow 


0.015 


9.5 


55 


0.11 



Table 3. The simulation: number of instabilities N, their size (in units of r jet ), their 
mean wave number k = 2irN/size, and the mean wavelength A = size/N for the multi-fast 

case. 













type 


N 


size 


k 


A 


Envelope 


3 


8.6 ± 0.75 


2.19 ± 0.19 


2.87 ± 0.25 


Core 


21 


10.5 ± 0.75 


12.57 ± 1.00 


0.50 ± 0.04 



Table 4. Wavelengths of the envelope X E , of the core Ac, and their ratio. 



method 




Ac 


ratio 


Stability analysis 


3.11 


0.53 


5.87 


Simulation 


2.87±0.25 


0.50±0.04 


5.74±0.96 



